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A  STUDY  OF  COMPRESSIBLE  TURBULENCE 


INTRODUCTION 

One  of  the  critical  aspects  of  the  aerodynamic  design  of  high  speed  aircraft  is  the  lack 
of  understanding  of  (compressible)  turbulence  at  these  speeds.  Turbulence  can  considerably 
affect  the  drag,  heat  transfer,  and  engine  performance  of  the  aircraft,  and,  consequently,  it  is 
imperative  that  the  understanding  of  turbulence  at  hypersonic  speeds  be  improved.  It  is 
proposed  that  as  a  first  step  the  differences  between  compressible  turbulence  and  the  much 
better  understood  incompressible  turbulence  be  better  identified.  In  Phase  I,  these 
differences  are  investigated  using  data  from  numerical  simulations  of  turbulence. 

Among  the  many  factors  that  influence  "hypersonic  turbulence”  are  compressibility,  the 
appearance  of  fluctuating  shock  waves,  real  gas  effects,  catalytic  wall  effects,  and  rarefied 
gas  effects.  The  influence  of  any  of  these  factors  on  turbulence  is  unknown,  and  it  would  be 
a  formidable  task  to  include  all  of  these  in  a  single  research  effort.  In  this  effort,  only 
compressibility  effects,  including  fluctuating  shock  waves,  are  considered.  This  factor  is 
common  to  all  hypersonic  turbulence  and  must  therefore  be  understood  before  substantial 
attention  is  devot^  to  the  other  factors. 

There  is  a  very  considerable  body  of  literature  on  turbulence  phenomena  at  low 
(incompressible)  speeds.^  These  investigations  of  low  speed  turbulence  are  either 
experimental,^  or  achieved  by  numerical  simulation.^  The  former  does  not  give  as  much 
detail  as  the  simulations,  but  does  give  information  at  realistic  Reynolds  numbers. 
Simulations  are  restricted  to  low  Reynolds  numbers  because  of  the  computer  limitations  and 
should  only  be  used  as  a  guide  to  the  nature  of  phenomena  at  higher  Reynolds  numbers.  At 
high  speeds  there  have  been  experimental  efforts^'^  to  obtain  data  on  hypersonic  turbulence, 
but  this  has  not  progressed  to  the  stage  where  underlying  phenomena  are  understood. 
There  is  a  general  belief^  that  compressibility  effects  in  turbulence  will  occur  only  when  the 
Mach  number  of  turbulence  fluctuations  achieves  magnitudes  of  0.5  or  higher.  The 
freestream  Mach  (M^  number  where  this  may  occur  depends  on  the  type  of  flow;  for  a 
boundary  layer,  M„  is  assumed  to  be  about  5,  whereas  for  mixing  layers,  it  is  about  1.5. 
Hypersonic  flight  vehicle  configurations  of  interest  may  fly  at  Mach  numbers  from  5  to  25 
depending  on  the  mission,  and  since  these  are  above  the  bounds  where  compressibility 
effects  become  important,  it  is  critical  that  the  nature  of  turbulence  at  these  speeds  be 
understood. 

Some  aspects  of  high  speed  turbulence  are  reasonably  well  understood.  For  example,  a 
concentrated  effort  in  compressible  shear  layers  over  the  past  several  years  has  identified 
mechanisms  which  cause  the  onset  of  reduced  turbulent  mixing  as  the  convective  Mach 
number  increases.  In  attached  boundary  layer  flows,  the  Van  Driest  II  theory  which 
accounts  for  near  wall  heating  performs  quite  well  at  freestream  Mach  numbers  up  to  at 
least  M„  -  5.0. 

However,  some  turbulent  flows  at  very  high  Mach  numbers  have  been  observed  to 
behave  differently  than  at  merely  high  Mach  numbers.  For  example,  the  laminar-to- 
turbulent  transition  in  boundary  layers  can  cover  extended  regions  of  flow,  rather  than 
being  a  brief  intermediate  stage  at  lower  speeds.  Furthermore,  earlier  studies^ ^  on  a 
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boundary  layer  indicate  that  pressure  fluctuations  dominate  the  turbulence  at  a  Mach 
number  of  9.4,  while  vorticity  fluctuations  dominate  low  speed  turbulence.  Hence,  there  are 
indications  that  turbulence  at  very  high  Mach  numbers  may  be  fundamentally  different  than 
at  low  or  moderate  Mach  numbers.  The  principal  distinguishing  characteristic  of  high  Mach 
number  flows  is  that  the  kinetic  energy  of  the  (free  stream)  flow  is  very  much  larger  than  the 
internal  energy  of  the  gas.  At  a  Mach  number  of  10,  for  example,  internal  energy  is  only  5% 
of  the  kinetic  energy.  High  speed  flows  are  affected  by  strong  near- wall  heating  resulting  in 
reduced  Reynolds  .  imbers  and  expanded  viscous  regions,  significant  effects  due  to  energy 
transferring  between  kinetic  energy  and  gas  compression,  and  shock  waves  which  are  often 
present. 

Because  compressibility  effects  appear  there  at  relatively  low  Mach  numbers,  one  of  the 
most  commonly  studied  problems^'°'^'^®  in  high  speed  turbulence  is  that  of  the  free  shear 
layer.  Here  two  streams  at  different  Mach  numbers  mix,  with  the  dominant  parameter  being 
a  convective  Mach  number,  M^.,  that  is  essentially  a  density  weighted  difference  of  the  Mach 
numbers  of  the  two  streams.  The  general  conclusion  of  these  studies  is  that  as  M^  increases, 
mixing  and,  by  implication,  turbulence,  decreases.  These  investigations  of  mixing  layers 
were  performed  initially  by  means  of  experiment,®  but  the  major  effects  of  compressibility 
have  also  been  delineated  by  numerical  simulation^  and  analysis.^*^  All  of  the  investigations 
are  for  M^  <  2,  which  represents  the  high  supersonic  regime. 

The  supersonic  mixing  layer  is  thus  an  accessible  problem  that  displays  compressibility 
effects  on  turbulence  at  relatively  low  speeds.  It  is  easily  computed  and  amenable  to 
mathematical  analysis,  and  thus  an  investigation  of  the  differences  between  turbulence  in 
compressible  and  incompressible  layers  should  be  straightforward.  In  particular,  data  from 
numerical  simulations  of  compressible  layers  contain  the  effects  of  incipient  shock 
formation,  a  feature  that  is  generally  believed  to  radically  alter  the  structure  of  turbulence. 
For  example,  it  is  suggested  in  Ref  10  that  just  as  a  shock  is  about  to  form,  the  nominally 
two-dimensional  mixing  layer  becomes  three  dimensional. 

An  important  part  of  the  proposed  work  is  to  investigate  the  role  of  vorticity  in 
compressible  turbulence.  Much  of  incompressible  turbulence  phenomena  can  be  described 
in  terms  of  vorticity  dynamics.  As  shown  later,  the  vortical  structure  of  turbulence  at  high 
Mach  numbers  is  considerably  different  than  at  low  speeds.  Density  and  pressure 
fluctuations  also  become  major  factors  at  high  speeds.  In  some  respects,  this  conclusion  is 
also  corroborated  in  the  mixing  layer  study  of  Ref.  10  which  implies  that  at  higher  M^ 
vortical  effects  are  reduced  because  of  density  changes. 

This  work  is  concerned  with  understanding  the  differences  between  incompressible 
and  hypersonic  turbulence  using  simulations  of  mixing  layer  flow  at  as  high  a  convective 
Mach  number  as  possible  so  that  dominant  compressibility  effects  can  be  clearly  identified. 

Objectives  -  Compressible  mixing  layers  will  be  examined  at  as  high  a  convective 
Mach  number  as  possible  in  order  to  delineate  compressibility  effects  on  turbulence.  The 
technical  questions  that  will  be  addressed  are  as  follows: 

(a)  What  are  the  differences  between  incompressible  and  hypersonic 
turbulence,  especially  the  role  of  vorticity  in  the  latter? 

(b)  What  is  the  role  of  shock  waves  or  shocklets  in  hypersonic  turbulence? 
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(c)  Is  the  standard  description  of  turbulence  based  on  velocity  fluctuations 
suitable  at  hypersonic  speeds  or  is  there  an  alternative  representation? 

Because  of  the  limited  time  and  costs  of  a  Phase  I  effort,  it  is  expected  that  only  a 
gross  understanding  of  these  topics  can  be  obtained.  This  should  provide  sufficient 
information  to  identify  areas  for  a  more  definitive  study  in  Phase  11. 

Principal  Results  -  This  work  involves  theoretical  analysis  of  turbulence  in  high 
speed  flow  and  large  eddy  simulations  of  the  planar  mixing  layer.  Some  of  the  key 
findings  are  as  follows.  One  analysis  indicates  that,  for  turbulence  to  exist  in  the  limit 
of  the  Mach  number  approaching  infinity,  turbulence  must  be  dominated  by 
streamwise  vortices.  An  approximate  model  based  on  the  idea  of  swept  vortices  at 
high  mach  number  makes  predictions  about  the  sweep  angle  of  these  vortices  from 
spanwise  at  low  speeds  to  streamwise  at  high  speeds,  and  about  the  reduced 
spreading  rate  at  high  speeds.  Simulations  of  planar  shear  lavers,  started  from 
random  initial  disturbances,  validate  the  structural  predictions  of  these  theories.  The 
simulations  predict  the  dominance  of  spanwise  vortical  structures  in  the  low  speed 
layer,  and  that  these  structures  are  swept  in  the  streamwise  direction  at  high  speeds. 
The  normalized  mixing  rate  is  predicted  to  decrease  significantly  at  high  speeds. 
Hence,  theory  and  simulation  make  consistent  predictions  about  the  structure  of 
turbulence.  There  is,  however,  a  significant  element  of  ambiguity  in  the  simulation 
results:  calculations  initiated  with  a  deterministic  disturbance  derived  from  stability 
theory  yielded  a  significantly  different  turbulent  structure. 

A  preliminary  analysis  has  been  performed  on  the  simulation  results,  and  some 
interesting  characteristics  about  the  turbulence  have  been  identified.  For  one, 
shocklets  are  generally  rare  and  weak,  even  at  a  convective  Mach  number  of  M^  =  2.5. 
The  turbulence  kinetic  energy  is  dominated  by  streamwise  fluctuations,  while  the 
other  two  components  of  energy  are  smaller  by  roughly  an  order  of  magnitude.  The 
pressure-velocity  correlations  promote  weak  transfer  of  energy  from  the  streamwise 
fluctuations  to  the  other  components  of  energy,  and  they  strongly  suppress  the  shear 
stress.  These  findings  are  compatible  with  the  unstead'-  pressure  and  vorticity  fields 
observed  in  the  simulations. 

Hence,  t.^is  rather  brief  Phase  1  effort  has  demonstrated  theory  and  simulations 
which  make  consistent  statements  about  high  speed  turbulence,  and  which  agree  with 
experimental  observations.  The  statistical  data  analyzed  so  far  indicate  dramatic 
changes  in  the  statistics  between  low  and  high  speed  flows  and  that  this  change  mav 
be  caused  bv  the  altered  vortex  structure  at  high  speeds. 

THEORY 

The  analysis  performed  during  this  effort  was  intended  to  determine  if  broad 
statements  can  be  made  about  the  nature  of  turbulence  as  the  Mach  number  becomes 
very  large.  The  thermodynamics  laws  which  lead  to  compressibility  effects  in 
turbulence  are  nonlinear  and  difficult  to  analvze.  However,  these  laws  can  be 
approximated  by  more  tractable  forms  over  limited  ranges  to  provide  insights  into 
phenomena  which  may  occur  in  high  speed  turbulence. 
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Inviscid  Model  -  It  is  frequently  assumed  that  large  scale  turbulence  can  be 
represented  by  the  Euler  equations;  although  viscosity  is  necessary  to  initiate  and 
sustain  turbulence,  its  global  features  can  be  represented  approximately  in  an  inviscid 
model.  The  density,  p,  temperature,  T,  and  pressure,  p,  in  the  steady  Euler  equations 
may  be  related  to  the  velocity,  q,  and  the  entropy  change,  S,  by 
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where  the  subscript  "r"  denotes  a  reference  value,  perhaps  in  the  region  of  turbulence 
production,  y  is  the  ratio  of  specific  heats,  and  R  is  the  gas  constant.  The  velocity  q  is 
normalized  with  respect  to  a  reference  value,  q,.. 


Let  q’  be  a  turbulent  fluctuation  about  q^,  that  is 

q^  =  1  +  2[Uj.u'  +  v^v'  +  w^w']  +  u'^  -t-  v'-  +  w'^ 


(4) 


and  let 


P/Pr  =  1  +  P' 

T/T,  =  1  +  T' 

p/Pr  =  1  +  P' 

S/R  =  S7R 

If  the  reasonable  assumption  that  the  general  fluctuating  quantity  f'  satisfies 

1  f'  I  <  <  1 


(5) 


(6) 


is  made,  then,  in  order  for  p/p^  to  be  real 


(y-1) 
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(q^-1) 


<  1 


(7) 


or 


I  u^u' +  v^v' +  w^w' +  (u'^  +  v'^  +  w'^)/2 1  <  l/[(y-l)M^] 


(8) 
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If  Eq.  (6)  is  applied,  Eq.  (8)  becomes,  to  first  approximation, 

lu^u'  +  v^v'  +  Wj-v/  I  <  l/[(y-l)M^]  (9) 

If  the  flow  is  such  that  and  w^,  then  Eq.  (9)  indicates  that 

lu^u'l  -♦O  as  (10) 

which  shows  that  the  fluctuations  become  two-dimensional  in  the  plane  normal  to  the 
streamwise  direction. 

In  a  similar  fashion,  the  relative  magnitude  of  the  nondimensional  fluctuations  in 
density,  pressure,  and  temperature  compared  to  the  velocity  can  be  determined.  Expanding 


Eq.  (1)  in  a  series  yields 

p'  =  -S7R-(qM)M2/2  +  0{[0.5(y-l)M~a-q2)]2}  (11) 

Then,  because  of  Eq.  (7),  a  first  approximation  to  Eq.  (11)  is 

p'  =  -S7R-(qM)M2/2  (12) 

From  Crocco's  equation,  entropy  is  related  to  vorticity  O  by 

TV(S/R)  =  (qxn)(yM2)  (13) 

in  which  q  and  T  are  the  velocity  vector  and  temperature.  By  a  rotational  perturbation 
about  an  irrotational  flow,  it  can  be  observed  that  to  first  approximation, 

Vx(S7R)  =  (qrXn')(yM2)  (14) 

Using  Eqs.  (12)  and  (14),  it  can  be  shovsm  that 

Vp'  -  -V(S7R)  -  m2  V(q2-l)/2  =  -m2  [yqxn'  +  V(q2-l)/2]  (15) 

Since  p'  is  bounded,  it  can  be  inferred  from  Eq.  (15)  that  as  M^  -» o°, 

y  qxfT  +  V(q2-l)/2  =  0  (16) 

If  q^  is  dominated  by  u^,  then  to  first  approximation,  Eq.  (16)  gives 

Uf  8u7 dx  =  0 

Ur  3u73y  =  yu^n^  (17) 

Ur  8u73z  =  -yUrfi^ 

Consider  now  the  vorticity  transport  equation 
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(18) 


L(n )  =  +  (q  •  V)n  -  (n-  V)q  +  n(V-  q)  =  -Vx(Vp/p)/ 

The  term  on  the  right  can  be  evaluated  by  combining  Eqs.  (1)  and  (3)  to  give 

-Vx(Vp/p)/  yM^  =  (y-1)  V(S/R)  x  Vq^  /  2  (19) 

or,  by  using  Eq.  (13), 

-Vx(Vp/p)/  yM^  =  (y-1)  (yM^/T)  (qxfl)  x  Vq2  /  2  (20) 

Using  Eq.  (20)  in  Eq.  (18)  and  assuming  that  the  normalized  vorticity  and  velocity  are 
bounded,  then  it  may  be  seen  that  as  -» 

(qxn)xVq2-»0  (21) 


If  q  is  dominated  by  u^,  then  Eq.  (21)  gives  the  following  possibilities: 


n  =0 

(22a) 

Urn2  =  Urn3  =  0 

(22b) 

u^  CI2  dq^/dy  -  u^  ^3  dq’/ dz  =  0 

(22c) 

Equation  (22a)  states  there  is  zero  vorticity,  and  Eq.  (22b)  indicates  that  streamwise  vortices 
dominate,  that  is,  a  Beltrami  flow  exists.  If  the  problem  is  a  geometrically  two-dimensional 

flow,  such  that 

1 dq-/ 3z 1  «  1 dq^/dy 1 

(23) 

then  Eq.  (22c)  indicates  that 

u^n2  dy  =  0 

(24) 

or 

02  =  0 

(25) 

This  indicates  that  vortices  may  either  be  streamwise  or  spanwise. 
combined,  it  may  be  seen  that 

If  Eqs.  (10)  and  (17)  are 

fb  =  03  =  0 

(26) 

leaving  the  only  possibility  that  the  vortices  are  predominantly  in  the  streamwise  direction. 

Physically,  the  analysis  indicates  that  once  vorticity  is  introduced  into  the  flow,  its 
distribution  is  dominated  by  the  entropy  relation  which  restricts  the  type  of  vorticity 
allowed.  To  a  first  approximation,  as  the  Mach  number  tends  to  infinity  in  a  typical  shear 
flow,  the  vortices  are  aligned  with  the  mean  streamwise  direction,  a  direction  which  causes 
no  entropy  production. 
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Simple  Vortex  Model  -  From  the  above  analysis,  it  appears  that  the  vortical  structures 
of  turbulence  will  tend  to  be  oriented  in  a  streamwise  direction  as  M  -» «>.  As  will  be  seen  in 
the  simulation  results  below,  the  structures  at  intermediate  Mach  numbers  are  highly  swept 
towards  the  streamwise  direction.  It  is  possible,  therefore,  that  the  sweep  angle  of  the 
structures  increases  with  increasing  Mach  number,  from  zero  sweep  at  incompressible 
speeds,  to  90°  as  M^.  -» 0°.  If  this  is  so,  how  does  this  structural  change  depend  on  the  Mach 
number? 

A  simple  theory  for  the  sweep  angle  of  the  dominant  structures  can  be  constructed  by 
postulating  that  the  vortical  structures  are  related  to  the  spanwise  ones  in  low  speed  shear 
layers,  but  that  they  respond  only  to  the  component  of  the  flow  which  is  normal  to  the 
vortices.  For  a  vortex  swept  at  an  angle  p  with  respect  to  the  spanwise  direction,  the 
"effective "  normal  velocity  difference  across  the  layer  and  convective  Mach  number  are 


AUgff  =  AU  cos(^) 

(27a) 

^c.cff  =  OOS(^) 

(27b) 

This  model  requires  the  empirically  derived  turbulent  growth  rate  for  two-dimensional 
structures 


<T2d(Mc)  (28) 

Take  (T2q  to  be  normalized  by  the  incompressible  growth  rate,  so  that  0’2d((’^  ~  From  two- 
dimensional  simulations,  cr2D  approaches  zero  as  in  Fig.  1  at  about  M^  =  1.  For  simplicity, 
this  function  is  approximated  as  (r2D  -  ^  ~  ^c-  This  approximation  may  be  least  valid  as 
Mj.  -♦  1,  but  this  is  largely  irrelevant,  it  turns  out,  to  the  predictions  made  by  the  model.  This 
analysis  implicitly  assumes  that  the  growth  rate  of  the  turbulent  layer  is  determined  by  the 
maximum  possible  growth  rate  over  all  sweep  angles,  which  seems  reasonable  but  is  not 
proven. 

Low  speed  shear  layer  experiments  have  demonstrated  that  the  layer's  growth  rate 
depends  upon  AU/U  in  which  U  is  the  average  of  the  velocities  on  either  side  of  the  shear 
layer^  By  the  above  postulate,  this  form  applies  for  high  speed  layers  with  AU^ff  replacing 
AU.  U  is  unchanged  bv  the  sweep  angle.  Hence,  the  velocity  difference  which  provides  the 
energy  for  the  growth  of  the  swept  structures  is  reduced  bv  a  factor  of  cos(^),  while  the 
convective  Mach  number  which  suppresses  growth  is  also  reduced.  The  growth  rate  of  the 
shear  laver  is  that  of  the  structures  with  no  restraint  on  sweep  angle,  which  is  predicted  to 
be 


^.id(  =  cos(^)  (TiqI  M^.  cos(^)  )  (29) 

IS  given  in  Fig.  2  for  several  angles  p.  At  low  Mach  numbers  the  maximum  growth  rate 
is  achieved  with  unswept  structures,  and  (T^plMj.,^)  is  unchanged  from  (T2D-  However,  at 
Mj,  -  0.6  the  swept  structures  have  a  more  rapid  growth  rate  than  the  unswept  ones.  As 
Mj,  ->  1  and  (T-i  D  (T  the  maximum  growth  rate  is  associated  with  increasingly  swept 
structures.  The  envelope  of  the  curves  for  a3Q  in  Fig.  2  is  similar  in  character  to 
experimental  measured  dependence  of  spreading  rate  on  the  Vlach  number. 


It  is  also  possible  to  determine  the  sweep  angle  at  which  the  growth  rate  is  a  maximum, 
as  a  function  of  convective  Mach  number,  plotted  in  Fig.  .3.  It  shows 

Pa,max  ~  ^  up  to  -  0.6,  at  which  point  increases  rapidly  with  M^,,  and 

asymptotes  to  as  It  is  also  worth  noting  that,  to  preserve  zero  net 

streamwise  angular  momentum,  there  must  be  vortices  with  opposite  signs  of  sweep  and 
streamwise  vorticity. 

It  is  worth  reiterating  that  this  is  a  very  simple  theory,  in  which  "turbulence"  is  highly 
idealized  and  even  then  treated  approximately.  Small  errors  are  to  be  expected;  it  is  general 
consistency  with  experiment,  simulation,  and  other  theory  that  is  sought. 

Other  theories  for  vortex  sweeping  at  high  convective  Mach  numbers  also  exist.  The 
theory  of  Nixon  and  Keefe^*^  postulates  that  the  vortex  structure  bifurcates  and  sweeps  to 
avoid  the  formation  of  a  shock  wave.  Typically,  shock  formation  occurs  at  M^  -  0.7  in  two- 
dimensional  simulations,  so  that  swept  structures  occur  for  M^  >  0.7.  Hence,  the  Nixon  - 
Keefe^^  theory  and  the  present  one  make  comparable  predictions  about  the  onset  of  swept 
structures,  but  the  underlying  mechanisms  appear  to  be  very  different.  Shocks  are  related 
to,  but  in  no  sense  directly  linked  to,  the  compressibility  effects  which  suppress  (T2d  with 
increasing  M^..  These  differences  may  become  critical.  If  a  good  model  is  to  be  constructed 
for  high  speed  turbulence,  it  is  important  for  the  model  to  be  based  on  the  correct 
underlying  physics.  Hence,  it  may  be  important  in  the  near  future  to  determine  which  of 
these  and  other  theories  gives  the  most  accurate  representation  of  the  phenomena  which 
control  the  turbulent  structures. 

Some  of  the  strongest  experimental  support  for  the  dominance  of  swept  structures  in 
high  speed  layers  comes  from  the  field  of  high  speed  jet  aeroacoustics.^’  The  structure  of 
the  dominant  turbulent  structures  in  the  near  field  of  a  round  jet,  whether  spanwise  or 
swept  vortices,  can  be  deduced  from  the  acoustic  pattern.  Axisymmetric  acoustic  modes  are 
caused  bv  axisymmetric  structures  which  are  spanwise  in  the  initial  jet  shear  layer,  while 
helical  modes  indicate  swept  structures.  Experiments  reveal  that  axisymmetric  modes 
dominate  the  near  field  up  to  jet  Mach  numbers  of  about  Mj  =  1.2,  which  corresponds  to 
Mj.  =  0.6.  Higher  speed  jets  contain  a  variety  of  helical  modes,  depending  upon  Mj.  While 
the  axisymmetric  shear  layer  is  different  than  the  plane  one  treated  in  the  above  theories, 
the  underlying  physical  mechanisms  are  similar.  Hence,  there  is  strong  experimental 
support  for  a  transition  to  swept  vortical  structures  at  M^  >  0.6  -  0.7. 

The  above  two  theories  described  in  Eqs.  (1)  -  (26)  and  in  (27)  -  (29)  differ  vastly  in  the 
assumptions  employed  and  the  depth  of  the  analysis.  However,  both  indicate  that  the 
vortices  in  a  turbulent  shear  layer  become  aligned  with  the  streamwise  direction  as  the 
Mach  number  becomes  very  large. 

TURBULENCE  SIMULATIONS 

At  this  point  the  focus  of  the  wurk  shifts  to  simulations  of  turbulence  which  will  lend 
strong  support  to  the  above  theoretical  findings  and  will  provide  further  insight  into 
turbulence. 
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Navier-SiOkes  Solver  -  TMRC  -  The  code  employed  to  perform  the  large  eddy 
simulations  (LES)  of  the  shear  laver  is  called  TMRC  (for  turbulence  modeling  research  code). 
A  brief  description  of  the  code  and  the  details  of  the  simulations  are  provided  here. 

The  equations  solved  by  TMRC  are  the  full  compressible  Navier-Stokes  equations  for  a 
perfect  gas  (y  =  1.4)  and,  if  appropriate,  equations  for  an  eddy-viscosity  or  Reynolds-stress 
turbulence  model.  The  Sutherland  formula  is  used  to  determine  the  laminar  viscosity. 
TMRC  uses  true  fourth  order  central  finite  volume  discretization  for  the  convective  terms 
and  second  order  central  differencing  for  the  viscous  fluxes.  Artificial  dissipation  is 
employed  in  high-Reynolds  number  and  in  high  Mach  number  flows  to  control  ''odd-even' 
oscillations  that  are  not  damped  by  viscosity,  especially  those  which  aj  'ear  at  shocks  or 
during  vortex  stretching  in  turbulent  flows.  Both  physical  phenomen,  'ause  a  rapid 
cascade  of  energy  to  high  wave  numbers  and  can  cause  the  solution  to  diverge.  A 
combination  of  linear  fourth  order  and  a  nonlinear  second/fourth  order  dissipation  is  used. 
In  the  nonlinear  dissipation,  second  order  viscous-like  dissipation  is  employed,  and  the 
coefficient  of  dissipation  depends  on  fourth  derivatives.  Fourth  derivatives  depend  most 
strongly  on  the  high  wave  number  content  of  the  solution;  hence,  the  nonlinear  dissipation 
behaves  like  a  dynamical  sub-grid-scale  model,  providing  diffusive  transport  only  if  there  is 
high  wave  number  content  in  the  solution.  An  explicit  five  stage  Rung^'  -  Kutta  procedure  is 
used  to  integrate  the  equations  in  time. 

The  effects  of  sub-grid-scale  (SGS)  turbulence  are  modeled  with  a  modified 
Smagorinsky  eddy-viscosity  model.  This  is  generally  to  be  regarded  as  a  poor  model  which 
has  several  shortcomings  for  the  present  flow.  However,  the  model  directly  affects  only  the 
smallest  scale  of  turbulence,  and  it  should  not  significantly  alter  the  dominant  features  in 
the  simulation  results. 

Boundary  Conditions  -  The  far-field  boundary  conditions,  which  must  permit 
acoustic  radiation  and,  in  one  case,  vortex  outflow,  are  treated  with  Thompson's*^  non¬ 
reflecting  boundary  condition  procedure.  These  boundary  conditions  are  reasonable  for 
acoustic  radiation,  but  very  poor  for  the  outflow  of  vortices.  Most  simulations  are  of  shear 
layers  which  are  periodic  in  both  directions  in  the  plane  of  the  shear  layer,  and  evolving  in 
time.  The  assumption  of  temporally  evolving  shear  layers  was  employed  to  provide  good 
resolution  of  turbulence,  within  the  limited  scope  of  the  simulations  that  can  be  afforded 
during  a  Phase  1  effort.  Given  the  assumption  of  periodicity,  these  boundaries  are  then 
treated  exactly,  and  unsteady  vorticity  and  acoustic  signals  pass  through  these  boundaries 
properly. 

It  is  important  to  insure  that  the  assumption  of  periodicity  does  not  compromise  the 
accuracy  of  the  simulation,  and  the  issues  of  major  concern  are  addressed  here.  One 
significant  concern  is  that  the  periodic  boundary  conditions  may  impose  significant 
constraints  on  the  structure  of  turbulence  which  would  render  tlie  simulated  iw  very 
different  from  the  natural  one.  The  computational  domain  used  here  has  an  a.  pect  ratio  in 
the  plane  of  the  shear  layer  of  where  x  is  the  streamwise 

direction,  z  is  spanwise,  and  y  is  normal  to  the  plane  of  the  shear  layer.  A  whole  number  of 
turbulent  structures  must  fit  into  this  domain  because  of  the  periodic  boundary  conditions, 
and,  typically,  a  small  number  of  structures  exist,  1,  2,  or  3.  These  conditions  constrain  the 
swept  structures  in  the  high  speed  shear  layer  to  certain  discrete  sweep  angles.  The  fewer 
the  number  of  structures  present,  the  larger  the  jumps  between  the  permitted  sweep  angles. 
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A  straightforward  means  of  reducing  this  problem  is  to  increase  the  size  of  the  domain  so 
that  a  large  number  of  structures  exist  in  the  domain. 

Simulations  were  perfo:  med  on  two  sizes  of  computational  domains.  One  was  50% 
larger  than  the  other  in  all  spatial  directions,  while  the  grid  spacing  and  the  initial  shear 
layer  thickness  were  unchanged.  The  principal  difference  between  the  results  was  that  the 
statistics  from  the  small-grid  solution  were  less  smooth.  Qualitatively,  the  dominant 
features  such  as  the  structures  and  the  growth  rate  were  similar.  This  result  suggests  that 
these  types  of  difficulties  associated  with  periodicity  are  under  control,  though  not 
necessarily  negligible.  All  simulations  reported  herein  were  run  on  the  larger  domain. 

Another  issue  with  the  periodic  assumption  is  that  of  entropy  and  vorticity  generation 
due  to  shock  waves  in  the  non-turbulent  flow  adjacent  to  the  shear  layer.  Turbulent  shear 
flows  are  clearly  not  periodic  in  the  streamwise  direction;  entropy  is  generated  if  shocks  are 
present,  and  turbulent  mixing  and  dissipation  generates  entropy.  However,  computational 
studies  contrasting  spatially  evolving  and  temporally  evolving  (streamwise  periodic)  low 
speed  shear  layers  exhibit  only  very  minor  differences.  A  comparison  of  temporally  and 
spatially  evolving  high  speed  layers  was  performed  in  support  of  this  work.  Because  of  the 
specific  initial  disturbances  used  to  excite  the  layer,  strong  shocks  formed.  However,  the 
vortical  structures  in  the  spatially  and  temporally  evolving  calculations  were  still  very 
similar.  Hence,  the  periodic  assumption  seems  good  enough  to  study  many  aspects  of  these 
flows.  In  fact,  shock  waves  were  seldom  observed  in  the  majority  of  the  simulations  which 
were  started  with  random  initial  perturbations. 

Code  Validation  -  The  code  TMRC  has  been  used  in  a  very  wide  variety  of  flow 
simulations,  and  its  ability  to  accurately  solve  the  Navier-Stokes  equations  has  been 
established.  Some  examples  of  previous  calculations  are  given  here.  A  pair  of  adjacent 
vortices  interact  to  spin  around  one  another,  and  this  behavior  can  be  described  analytically 
in  the  limit  of  negligible  viscous  effects.  Calculations  of  these  flows  were  performed,^'^  and 
the  computed  solution  converges  to  the  analytical  one  as  the  grid  is  refined  and  the  viscosity 
goes  to  zero.  Fourth  order  accuracy  was  demonstrated.  The  code  has  been  used  to  perform 
simulations  of  the  axisymmetric  mean  flow  and  turbulence  structure  of  a  round  jet.  In  the 
near  field  where  the  assumption  of  axisymmetric  turbulence  is  reasonable,  the  predicted 
shock  structure  and  mean  velocity  profiles  agreed  well  with  experimental  measurements.^'* 
The  code  has  also  been  used  to  perform  large  eddy  simulations  of  stratified  turbulence,*^ 
and  simulations  of  a  variety  of  impinging  jet  phenomena.*^'* '  In  all  cases  the  code  has  given 
results  that  agree  well  with  analysis  or  experiment  when  limiting  assumptions  are  satisfied. 
Hence,  there  is  a  high  level  of  confidence  in  the  solutions  produced  by  TMRC. 

TMRC  has  also  been  used  to  examine  the  enhanced  turbulence  in  a  jet  impinging  on  the 
ground.  An  intriguing  result  of  this  study  was  the  formation  of  fluctuating  shock  waves  (or 
shocklets)  due  to  the  turbulent  eddies  in  the  impingement  region.  The  simulations  correctly 
predicted  the  enhancement  of  turbulence  at  specific  frequencies  due  to  fluid/acoustic 
feedback  and  unsteady  shock  formation. 

Domain  and  Grid  -  The  grids  used  for  these  simulations  are  Cartesian.  For  the 
spatially  periodic  simulations,  the  grid  spacing  in  the  x-  and  z-directions,  which  are  in  the 
plane  of  the  shear  layer,  is  uniform  with  Ax/6q  =  Az/6g  =  0.5,  in  which  6o  is  the  initial  shear 
layer  thickness.  Uniform  grid  spacing  in  the  x-  and  z-directions  was  used  to  minimize  any 
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bias  in  the  formation  of  spanwise  or  streamwise  turbulent  structures.  In  the  y-direction, 
normal  to  the  plane  of  the  shear  layer,  the  grid  is  stretched,  with  a  grid  spacing  of 
Ay/6o  =  0.2  throughout  the  region  of  turbulence,  stretching  to  Ay/6o  =  1.5  at  edges  of  the 
domain.  The  extent  of  the  computational  domain  is 

-18<x/6o<18,  -19.2  <y/6o<  19.2,  and  -8.5  <  z/6o  <  8.5 

The  grid  used  for  the  majority  of  simulations  had  73  x  85  x  35  points  in  the  x-,  y-,  and  z- 
directions.  Statistics  and  structural  data  are  obtained  from  the  shear  layers  when  they  are 
typically  about  6/6q  =  4,  so  that  the  dimensions  of  the  domain  are  reasonably  large  relative 
to  the  shear  layer  thickness. 

FLOW  FIELDS 

Free  shear  layers  only  were  considered  in  this  work.  These  are  relatively  simple  flows, 
which  are  well  studied,  and  which  exhibit  strong  compressibility  effects.  The  approach 
taken  here  is  not  limited  to  free  shear  flows;  it  can  also  be  applied  to  temporally  or  spatially 
evolving  boundary  layers  and  other  flows. 

Temporally  Evolving  Flows  -  The  principal  characteristics  which  define  each  flow  are 
the  mean  velocity  and  thermodynamic  variables,  and  the  initial  disturbance  field  which  will 
grow  and  become  turbulence.  The  mean  initial  field  consists  of  an  approximately  linear 
variation  between  the  two  end  states  of  the  velocity.  In  the  one  simulation  which  involved 
mean  density  variation,  the  density  also  varied  linearly.  The  initial  pressure  was  uniform  in 
all  cases. 

Two  types  of  initial  disturbance  fields  were  considered.  In  aU  but  two  simulations,  the 
initial  velocity  field  was  computed  as  a  smoothed  random  field.  Initially  the  perturbation 
velocity  field  is  specified  randomly  with  a  maximum  magnitude  of  AU,  the  total  velocity 
difference  across  the  layer.  The  Cray  random  number  generator  RANF  is  used  to  obtain  the 
random  numbers.  The  initial  disturbances  were  confined  to  a  layer  20%  thicker  than  the 
mean  velocity  profile,  and  they  were  scaled  to  be  a  maximum  on  the  centerplane  and  to 
decay  smoothiy  to  zero  at  the  edge  of  the  region  of  disturbances. 

u'  =  AU  RANF  F(  I  y-yp  I  /(0.66o) )  (31) 

The  scaling  function  goes  smoothly  from  F{0)  =  1  to  F(l)  =  0.  A  different  random  number  is 
used  for  each  component,  u',  v',  and  w',  of  the  vector,  u'.  This  u'  field  is  then  added  to  the 
mean  velocity,  and  the  complete  velocity  field  is  smoothed  with  second  order  smoothing, 
e.g.. 


Un+l  .  yn  ^  ^^2^^  (32) 

in  which  0"  denotes  the  iteration  and  e  is  a  small  number,  until  the  velocity  fluctuations  are 
reduced  considerably.  The  energy  spectrum  resulting  from  this  smoothing  process  has  a 
maximum  at  low  wave  numbers  and  decreases  with  increasing  wave  number,  as  occurs  in 
physical  turbulence.  This  velocity  field  is  unlike  low-level  turbulence,  and  the  first  process 
to  occur  during  the  simulation  is  for  the  flow  to  relax  to  a  physically  plausible,  although 
very  random,  disturbance  field. 
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These  initial  conditions  contain  no  organized  structures,  and  any  structures  that  appear 
in  the  solutions  are  formed  as  a  result  of  the  dynamical  processes  that  govern  turbulence. 

The  other  type  of  initial  disturbance  was  used  in  the  spatially  evolving  calculation  and 
a  companion  temporal  simulation.  It  was  a  combination  of  three  modes  selected  using  the 
results  of  stability  analyses.  The  fundamental  two-dimensional  mode  for  the  initial  shear 
layer,  its  first  subharmonic,  and  a  pair  of  oblique  waves  were  combined  for  this  case.  This 
initial  condition  for  the  temporally  evolving  flow  is  the  appropriate  equivalent  of  the  inflow 
conditions  supplied  to  the  spatially  evolving  flow  described  below. 

Spatially  Evolving  Flow  -  The  simulation  of  an  M^.  =  1.,  spatially  evolving  shear  layer 
was  conducted  on  a  401  x  76  xl3  (x,y,z)  grid  spanning  a  spatial  domain  of  0  <  x/6o  <  500, 
-15  <  y/6o  <  15,  and  -3  <  z/6q  <  3.  The  value  of  M^.  =  1.  was  obtained  by  imposing  M|=3.1 
and  M2  =1.1,  with  constant  density  and  pressure  across  the  layer  at  the  end  of  the  splitter 
plate.  The  layer  was  excited  at  the  plate  trailing  edge  by  spatial  instability  modes  calculated 
from  the  stability  code  of  Sandham  and  Reynolds’® .  To  the  most  unstable  2D  mode  and  its 
subharmonic  were  added  a  pair  of  oblique  waves  (^  =  78°)  that  had  growth  rates  very  near 
the  2D  disturbance.  For  these  waves  the  growth  rate  was  about  50%  of  the  maximum 
achievable  at  ^  =  60  °  for  this  convective  Mach  number.  The  simulation  was  run  for  times 
sufficient  to  allow  initial  disturbances  at  the  splitter  plate  to  pass  completely  out  of  the 
domain. 

SIMULATION  RESULTS 

The  principal  results  are  from  shear  layers  with  a  Mach  number  difference  across  the 
layers  of  AU/c  =  0.4  and  AU/c  =  5.0  with  a  uniform  mean  density.  The  flows  correspond  to 
nominal  convective  Mach  numbers  of  M^  =  0.2  and  M^.  =  2.5,  respectively.  In  a  third  case  the 
velocity  difference  AU  =  8.  The  density  was  varied  to  mimic  the  effects  of  constant 
stagnation  temperature  across  the  layer,  with  fast  cold  flow  above  the  shear  layer  and  hot 
stagnant  fluid  below  it.  This  gives  a  nominal  value  of  M^.  =  1.48,  assuming  the  convective 
Mach  numbers  are  the  same  with  respect  to  upper  and  lower  freestreams.  The  final 
simulation  to  be  discussed  was  at  AU  =  8  and  constant  density  to  give  M^.  =  4. 

Shear  Layer  Growth  -  The  simulations  were  run  until  the  shear  layer  thickness  grows 
to  6  -  68q,  typically.  The  growth  is  nonlinear  while  6  <  36^,  but  is  roughly  linear  for 
46q  <  8  <  68^.  Figure  4  gives  the  shear  layer  thickness  as  a  function  of  time  for  three 
different  simulations.  The  linear  growth  rate  suggests  that  the  turbulence  has  reached  some 
asymptotic  state,  beyond  the  effects  of  the  initial  conditions  and  before  the  structures 
become  saturated  in  the  computational  domain.  The  results  which  are  examined  in  detail 
and  presented  later  are  taken  from  this  linear  growth  period. 

There  is  a  long  initial  latency  period  in  the  two  high  speed  flows,  during  which  time 
the  growth  rate  is  even  slower  than  its  asymptotic  value.  This  may  be  due  to  the  initial 
conditions.  However,  the  normalized  initial  velocity  perturbations  were  similar  for  all  three 
cases  represented  on  this  figure,  and  this  result  may  reflect  a  significant  physical 
characteristic  of  developing  shear  layers.  There  has  been  no  fu.rther  investigation  of  this 
issue  during  the  present  work. 
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The  asymptotic  growth  rates  of  these  shear  layer  are  given  in  Fig.  5.  The  growth  rates 
are  normalized  by  an  incompressible  rate  which  is  obtained  by  assuming  the  M^.  =  0.2 
simulation  to  be  essentially  incompressible.  There  are  few  simulation  data  points,  and  the 
simulation  at  M^.  =  1.48  had  a  strong  density  gradient  across  the  layer.  However,  the  trend 
observed  in  the  simulations  is  clearly  toward  reduced  turbulent  mixing  at  high  speeds. 

Dominant  Struchme  -  There  are  significant  differences  in  the  vortex  structures  at  low 
and  high  speeds.  A  side  view  of  the  spanwise  vorticity  for  the  M,.  =  0.2  case  is  given  in 
Fig.  6-a.  The  vorticity  has  largely  coalesced  into  three  principal  "clumps."  Figure  6-b  gives  a 
top  view  of  pressure  contours  which  identify  low  pressure  regions  on  a  few  planes  near  the 
middle  of  the  shear  layer.  The  low  pressure  is  associated  with  the  vortex  cores.  The  now 
classical  picture  of  a  flow  dominated  by  spanwise  vortices  is  evident,  but  a  considerable 
amount  of  randomness  in  this  structure  is  also  present.  Recall  that  this  turbulence  field  is 
emerging  from  a  purely  random  initial  disturbance,  and  that  the  mean  shear  layer  thickness 
has  grown  by  a  factor  of  -4  only.  The  vorticity  plot  for  this  flow  is  not  used  to  illustrate 
structure.  The  vorticity  is  sensitive  to  high  wave  number  fluctuations  which  tend  to  be 
random,  and  this  produces  a  very  confusing  plot.  The  pressure  provides  a  form  of 
integration  over  all  fluctuations  and  gives  a  better  picture  of  the  large  scale  structure.  The 
result  demonstrates  that  the  flow  selects  spanwise  vortical  structures. 

The  swept  structures  in  the  M^.  =  2.5  case  can  be  seen  in  the  vorticity  contour  plots 
given  in  Fig.  7.  The  top  view  in  Fig.  7-a  gives  the  streamwise  vorticity  on  several  planes 
through  the  layer.  The  contour  levels  were  selected  to  highlight  the  high-vorticity  vortex 
cores.  It  was  necessary  to  plot  the  vorticity  on  several  planes  because  the  vortices  are  also 
vertically  inclined  at  a  shallow  angle.  This  vertical  inclination  is  seen  in  a  side  view  which 
shows  the  streamwise  velocity  in  a  single  x-y  plane.  The  top  view  indicates  that  the  vortices 
are  highly  swept,  at  roughly  p  =  77°  from  the  spanwise  direction.  An  oblique  view  from 
above  the  layer,  in  Fig.  7-c,  gives  streamwise  vorticity  on  z-y  planes  which  are  uniformly 
spaced  in  the  streamwise  direction.  The  highly  swept  vortices  are  again  seen.  Figure  8  gives 
a  top  view  of  the  contours  of  the  thermodynamic  pressure.  The  contour  levels  are  selected  to 
highlight  low  pressure  which  is  associated  with  the  vortex  cores  seen  in  Fig.  7-a.  The 
minimum  contour  level  is  0.45  of  the  freestream  pressure  but  the  minimum  pressure  in 
the  simulation  was  about  0.3  P„. 

The  majority  of  the  pressure  fluctuations  in  this  simulation  are  closely  related  to  the 
vortex  structure.  This  may  be  of  considerable  importance  because,  as  will  be  discussed 
below,  the  suppression  of  turbulence  can  be  viewed  as  an  effect  of  pressure-velocity 
correlations.  If  these  pressure-velocity  correlations  depend  on  the  vortex  structure,  then 
knowledge  of  the  vortex  structure  will  be  required  for  statistical  turbulence  prediction 
methods. 

As  discussed  in  the  section  on  boundary  conditions,  the  structure  of  turbulence  is 
constrained  bv  assumption  of  spatial  periodicitv.  The  turbulence  stmcture  in  the  M^.  =  2.5 
simulation  contains  two  fundamental  periods  in  the  spanwise  and  one  in  the  streamwise. 
This  simulation  would  force  turbulence  with  a  natural  structure  of  for  example,  1.5  or  2.5 
periods  in  the  spanwise  direction  onto  the  solution  given  in  Fig.  7.  Hence,  the  preferred 
sweep  angle  could  be  in  the  range  -70°  to  -80°.  This  ambiguity  does  not  alter  the 
conclusion  that  highly  swept  vortices  dominate  the  high  speed  mixing  layer. 
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The  effective  convective  Mach,  which  is  the  Mach  number  normal  to  the  structures,  is 
given  by  =  M^  cosip).  At  Mj.  =  2.5  and  p  =  77°,  M^  gf£  =  0.56.  This  is  clearly  a  subsonic 
Mach  number,  but  transonic  flow  could  be  generated,  depending  on  the  nature  of  the 
structure.  This  simulation  result  is  consistent  with  the  theories  presented  above;  The  sweep 
angle  is  very  large,  even  at  this  moderate  Mach  number,  and  of  roughly  the  correct 
magnitude.  At  M^.  =  0.56,  these  turbulent  structures  exist  just  below  the  Mach  number  at 

which  a  (further)  change  in  the  angle  of  the  structure  is  predicted  by  Eq.  (32)  and  in  Fig.  3. 

The  vortex  structure  in  the  high  speed  case  is  well  defined  and  it  exhibits  considerably 
less  randomness  than  in  the  low  speed  case.  The  sub-grid-scale  viscous  effects  used  in  the 
simulations  scales  linearly  with  Mach  number,  which  would  provide  a  similar  effective 
turbulent  Reynolds  number  if  the  normalized  turbulence  were  the  same  in  both  flows.  The 
less  vigorous  normalized  turbulence  in  the  high  speed  case  causes  the  effective  turbulent 
Reynolds  number  to  be  lower,  which  restricts  the  range  of  scales  of  turbulence  and  makes  it 
less  chaotic.  This  is  a  deficiency  which  can  easily  be  remedied,  now  that  it  is  recognized.  In 
the  present  results,  it  should  principally  affect  the  small  turbulent  scales,  while  the  largest 
scales  should  be  much  less  affected.  The  limited  range  of  scales  will  affect  the  predicted 
levels  of  anisotropy  because  the  small  scales  tend  to  be  more  isotropic,  while  the  larger 
scales  are  less  isotropic.  Hence,  the  predicted  levels  of  anisotropy  may  be  artificially  larger 
than  in  a  high  Reynolds  number  shear  layer. 

A  very  different  turbulent  structure  was  seen  in  two  of  the  simulations.  In  the  spatially 
and  temporally  evolving  flows  initialized  with  three  modes  selected  from  instability  theory, 
the  vortical  structure  was  essentially  spanwise.  While  the  flow  was  at  a  lower  convective 
Mach  number,  M^  =  1.0,  it  is  in  the  range  where  the  theory  predicts  a  sweep  angle  of 
p  -  60° .  The  lack  of  sweep  in  this  simulation  suggests  a  strong  dependence  on  the  initial 
disturbance.  It  is  possible  that  the  flow  is  non-unique,  a  characteristic  which  could  be 
exploited  in  flow  control.  However,  this  result  clearly  illustrates  a  potential  hazard  of 
initiating  the  flow  with  highly  regulated  disturbance  modes. 

With  the  above  paragraph  as  a  caveat  on  the  veracity  of  the  simulation  results,  it  must 
further  be  noted  that  Papamoschou  and  Roshko®  have  speculated  that  the  dominant 
structures  in  their  experiments  were  spanwise.  Their  argument  was  that  structures  were 
identified  in  optical  records  which  integrated  across  the  span  of  the  layer,  and  that  spanwise 
structures  only  would  be  preserved  through  the  integration  process.  That  speculation  and 
the  present  findings  are  clearly  in  disagreement. 

Shocklets  -  Two-dimensional  simulations  of  compressible  mixing  layers  have 
demonstrated  the  presence  of  shocklets  which  are  associated  with  the  passage  of  vortical 
structures.  Transient  shocks  have  also  been  observed  in  three-dimensional  simulations  of 
homogeneous  turbulence.  Furthermore,  the  existence  of  shocklets  is  central  to  some 
theories  which  attribute  the  reduced  turbulent  growth  rates  of  high  speed  shear  layers  to  the 
extra  dissipation  caused  by  these  shocklets. 

As  the  vortical  structures  become  skewed,  the  effective  Mach  number  decreases  as 
cosip).  The  vortex  structure  observed  in  the  M^.  =  2.5  case  was  swept  so  that  M^  =  0.56, 
and  shocklets  were,  almost  universally,  not  seen.  At  only  one  instant  in  time  was  a  shock 
observed  in  the  simulations  of  this  and  other  high  speed  flows  started  with  the  random 
initial  conditions.  An  effort  was  made  to  track  that  shock,  but  it  disappeared  very  quickly. 
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The  principal  question  would  then  seem  to  be:  Under  what  conditions  do  shocklets 
constitute  a  significant  phenomenon  in  turbulence?  Shocklets  were  present  in  both  spatially 
and  temporally  evolving  simulations  performed  with  the  starting  conditions  prescribed 
from  stability  theory.  Clearly,  the  simulations  are  capable  of  representing  shocklets,  and  a 
flow  with  shocklets  can  evolve.  Of  the  present  results,  those  started  from  random  initial 
conditions  would  seem  to  be  more  representative  of  a  "natural"  state  of  turbulence. 
However,  this  may  not  be  so  if  the  flow  is  prone  to  have  a  lengthy  "memory"  of  the  initial 
disturbances.  Furthermore,  certain  types  of  additional  strains  such  as  lateral  divergence 
(dW/dz  >  0,  in  the  present  notation)  would  reduce  the  sweep  angle  of  the  vortices,  literally, 
by  stretching  them  into  a  more  spanwise  posture.  This  would  increase  their  convective 
Mach  number  and  the  likelihood  of  shocklet  formation. 

One  result  of  the  present  work  is  that  shocklets  are  not  required  to  reduce  the  growth 
rate  of  a  shear  layer  at  high  convective  Mach  number.  The  significant  growth  rate 
suppression  in  the  M^  =  1.48  and  2.5  cases  was  caused  by  a  change  in  the  dominant  vortical 
structure,  in  the  absence  of  shocklets. 

=  4  Case  -  A  major  goal  of  this  work  was  to  examine  the  results  of  simulations  at 
as  high  a  Mach  number  as  possible.  The  M^.  =  2.5  case  falls  in  the  high  end  of  the  supersonic 
range,  but  a  case  at  M^,  =  4  is  into  the  realm  of  hypersonic  flow.  Calculations  were 
performed  using  the  same  grid  and  procedure  for  generating  the  initial  conditions  as  the 
M^  =  2.5  case.  However,  the  initial  disturbances  decayed  and  the  flow  showed  little 
tendency  to  generate  a  "turbulent"  flow.  Two  more  calculations  of  this  case  were 
performed,  each  with  a  reduced  coefficient  in  the  sub-grid-scale  diffusion  model,  but  this 
did  not  change  the  outcome. 

The  inability  to  generate  a  growing  turbulent  layer  is  most  likely  caused  by  some  form 
of  computational  restriction  imposed  on  the  solution,  but  the  implications  for  flow  physics 
cannot  be  completely  disregarded.  It  may  be  that  the  amplified  structures  are  so  highly 
inclined  that  they  did  not  "fit"  into  the  computational  domain  properly.  Since  the  M^.  =  2.5 
case  had  only  one  structure  in  the  streamwise  direction  and  the  ratio  of 
streamwise/spanwise  dimension  is  expected  to  increase  with  Mj.,  it  is  possible  that  a  longer 
computational  domain  might  be  required  to  obtain  a  turbulent  flow.  (These  possibilities 
were  not  recognized  until  the  analysis  was  nearly  completed,  and  there  was  insufficient  time 
to  attempt  further  simulations.) 

It  is  unfortunate,  but  the  failure  to  achieve  a  turbulent  solution  in  this  case,  for  reasons 
which  may  be  due  to  limitations  of  the  numerical  simulations,  does  cast  some  doubt  on  the 
accuracy  of  the  simulation  at  M^  =  2.5. 

Turbulence  Statistics  -  The  statistics  of  these  turbulent  flows  provide  a  completely 
different  way  to  look  at  exactly  the  same  physical  phenomena.  Statistics  are  important  to 
consider  because  all  mean  flow  calculations  rely  upon  some  form  of  statistical  treatment  of 
turbulence. 

Turbulence  statistics  are  obtained  by  averaging  over  the  x-  and  z-directions,  for  which 
the  turbulence  is  periodic,  at  one  instant  in  time.  The  convergence  of  the  statistics  depends 
on  the  x-z  area  of  the  domain,  and  for  these  calculations  the  statistics  are  reasonably  well. 
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but  not  completely,  converged.  That  is,  a  simulation  on  a  larger  domain  or  a  different 
realization  of  the  same  simulation  would  yield  a  slightly  different  answer.  This  is  verified 
by  noting  that  the  normalized  statistics  are  slightly  different  at  different  times,  and  that 
some  statistics  which  should  be  exactly  zero  are  slightly  non-zero.  However,  the  important 
trends  that  are  identified  in  the  data  are  unaffected  by  these  small  uncertainties. 

In  general,  mass  weighted  or  Favre'  averaging  is  used,  in  which 

u  =  U  +  u",  U  =  <pu>/<p>  (33) 

Conventional  Reynolds  averaging  is  also  employed  occasionally.  Averages  are  denoted  by 
<>  brackets,  although  the  "overbar”  notation  is  also  used  elsewhere  in  this  report  (e.g., 

<g>  =  g)- 


In  general  the  statistical  results  are  normalized.  Spatial  dimensions  are  normalized  by 
the  shear  layer  thickness  6,  for  example,  q  =  y/6.  The  edges  are  taken  to  be  the  points  at 
which 


U(yedge)  =  U„,in  +  0.1AU  and  UCy^dge)  =  '  0-1  AU  (34) 

The  turbulent  stresses  and  other  higher  order  statistics  are  normalized  by  appropriate 
combinations  of  AU  and  the  mean  density,  so  that  effects  of  the  differing  mean  velocities  is 
removed  from  the  results. 

In  the  following  discussion,  the  results  at  =  0.2  and  =  2.5  are  contrasted.  The  low 
Mach  number  case  is  virtually  incompressible,  while  the  high  speed  case  is  very  much 
affected  by  compressibility.  The  density  ratio  across  both  shear  layers  is  unity. 

The  first  data  to  consider  are  the  mean  momenta  and  density  profiles,  not  normalized, 
given  in  Fig.  9.  The  low  speed  case  exhibits  a  smooth  hyperbolic  tangent-like  streamwise 
velocity  profile,  while  p  and  pV  are  roughly  constant,  as  expected.  The  high  speed  case 
displays  a  streamwise  velocity  profile  which  has  rather  sharp  corners,  in  contrast  to  the 
smooth  profile  at  the  lower  Mach  number.  This  difference  in  the  velocity  profiles  is 
somewhat  unexpected,  and  such  pronounced  sharp  corners  are  not  generally  observed  in  a 
wide  range  experimental  data.  However,  this  trend  is  consistent  some  experimental 
measurements.  The  sharp  corners  on  the  velocity  profile  cause  the  maximum  velocity 
gradient  dU/ayn,ax  t»e  less  than  in  a  shear  layer  of  similar  total  thickness,  but  at  a  low 
convective  Mach  number.  This  causes  the  vorticity  thickness  6^,  defined  by 
AU  =  (dU/  dy)^3^  and  normalized  on  its  incompressible  value,  to  be  greater  at  higher  M^.. 
It  has  been  observed  experimentally^'®  that  the  normalized  growth  rate  for  6^^  is  twice  as 
large  as  the  "visual”  thickness  or  pitot  thickness  growth  rate  at  high  M^.  Hence,  the  trend 
toward  a  more  nearly  linear  velocity  profile  is  not  inconsistent  with  experimental 
measurements. 

The  mean  density  is  about  40%  less  than  its  freestream  value  near  the  center  of  the 
shear  layer.  At  the  center  of  the  layer,  turbulence  accounts  for  a  significant  fraction  of  the 
total  normal  stresses,  turbulence  plus  thermodynamic  pressure.  The  turbulent  <pv”v"> 
normal  stress  displaces  some  of  the  thermodynamic  pressure  and  the  density  decreases  with 
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the  pressure.  The  instantaneous  thermodynamic  pressure  in  the  cores  of  vortices  in  the 
shear  layer  is  observed  be  as  low  as  0.3  P„  in  the  M^.  =  2.5  case. 

The  components  of  the  normalized  turbulent  stress  tensor  <pUiUj’>  are  given  in  Fig.  10. 
In  the  low  speed  case,  the  maximum  normal  stresses  are  in  the  range  0.03  to  0.04  and  the 
maximum  magnitude  of  the  principal  shear  stress  is  about  -0.02.  These  stresses  are 
somewhat  higher  than  typically  observed  in  experiments.^'^  The  normal  stresses  are  high  by 
perhaps  20%  to  40%  while  the  shear  stress  is  roughly  60%  above  the  upper  end  of  the 
experimental  range.  The  cause  of  this  discrepancy  is  unknown.  It  may  be  due  to  the  effects 
of  the  limited  range  of  scales  in  the  simulation,  a  transitional  effect,  or  the  strength  of  the 
initial  disturbance.  This  is  an  issue  which  deserves  further  scrutiny.  The  <pu"w">  and 
<pv  "w  ">  shear  stresses  are  close  to  zero,  as  they  should  be. 

The  stresses  in  the  high  speed  case,  given  in  Fig.  10,  are  very  different  than  in  the  low 
speed  layer.  The  maximum  streamwise  normal  stress  <pu''u”>  is  close  to  the  value  in  the 
low  speed  case.  However,  the  other  normal  stresses  are  approximately  an  order  of 
magnitude  less  than  in  the  low  speed  layer.  Hence,  there  is  a  very  unequal  distribution  of 
turbulence  energy  among  the  normal  stresses,  with  the  majority  of  the  energy  in  the 
streamwise  component.  The  shear  stress  is  also  significantly  less  than  in  the  low  speed  case, 
which  is  of  course  necessary  to  be  consistent  with  the  reduced  growth  rate. 

The  accuracy  of  the  normal  stress  results  in  the  =  2.5  case  is  a  point  of  concern.  This 
level  of  anisotropy  is  virtually  unknown  in  any  low  speed  turbulence.  Only  in  transitional 
flows  is  turbulence  essentially  two-dimensional  and,  hence,  highly  anisotropic.  As  noted 
previously,  the  SGS  transport  does  scale  with  the  Mach  number,  which  will  cause  the 
Mj.  -  2.5  case  to  be  more  laminar-like,  and  which  may  artificially  increase  the  level  of 
anisotropy.  However,  the  global  features  of  the  =  2.5  shear  layer  are  consistent  with 
theory  and  experiment.  It  is  inappropriate  to  reject  one  result  from  this  simulation  while  so 
many  other  aspects  are  reasonably  accurate.  Hence,  the  trend  toward  higher  levels  of 
anisotropv  and  high  convective  Mach  numbers  should  be  believed  unless  disproven. 

One  of  the  goals  of  this  work  is  to  determine  the  origins  of  the  reduced  shear  stress,  in 
terms  of  turbulent  statistics.  One  possibility  is  that  of  a  "purely  compressible"  mechanism. 
In  a  "purely  compressible"  mechanism  the  instantaneous  density  fluctuations  due  to  fluid 
dynamic  compression  would  correlate  with  velocity  fluctuations  in  a  manner  that  would 
reduce  the  shear  stress.  This  possibility  can  be  assessed  by  decomposing  the  full  stress 

<pu'j’Uj>  (35) 

into  its  incompressible  part 

<p><u5u'>  (36) 

and  a  part  which  is  purely  due  to  compressibility  effects,  written  for  convenience  as 


<pu'i'Uj>  -  <p><u'u'>  (37) 

If  there  is  a  high  correlation  between  the  blobs  of  fluid  which  carry  the  fluctuating  second 
moment  of  velocity  uju-  and  the  perturbation  density,  then  the  purely  compressible  part. 
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Eq.  (37),  will  be  significant.  Figure  11  gives  this  decomposition  for  <pu"u">  and  <pu"v">.  In 
the  M(,  =0.2  case,  the  compressible  part  is  negligible,  as  expected.  At  M^.  =  2.5  the 
compressible  part  is  not  zero,  but  it  is  roughly  an  order  of  magnitude  less  than  the  total 
stress.  The  decomposition  for  the  shear  stress  <pu  "v">  at  M^.  =  2.5  also  shows  the 
compressible  part  to  be  small.  The  compressible  part  generally  has  the  opposite  sign  of  the 
stress  itself,  and  hence  contributes  to  a  reduction  in  the  stress.  For  this  case  in  which  there  is 
no  density  stratification,  the  total  stresses  <pu'j’Uj'>  are  not  very  different  from  the 
incompressible  part  of  the  stresses  <p><U'Uj>.  This  result  suggests  that  purely  compressible 
mechanisms  which  depend  upon  significant  correlations  between  the  instantaneous  velocity 
fluctuations,  u'uj,  and  the  perturbation  density  may  not  be  very  important. 

Anisotropy  of  the  normal  stresses  is,  however,  very  important.  The  principal  source  of 
the  shear  stress  is  its  production  term, 

=  -<pv"v">  dU/dy  (38) 

for  this  simple  shear  layer.  The  reduced  normal  stress  <pv”v">  almost  guarantees  that  shear 
stress  will  be  small  as  well,  unless  drastic  reductions  in  the  mechanism  which  destroy 
<pu  "v”>  also  occur.  Hence,  the  origins  of  the  anisotropic  normal  stresses  may  be,  from  the 
statistical  point  of  view,  a  important  contributor  to  compressibility  effects. 

Pressure- Velocity  Terms  -  Anisotropy  is  regulated  by  fluctuating  pressure-velocity 
interactions.  In  simple  shear  flows,  the  turbulence  energy  is  formed  as  streamwise 
fluctuations.  This  energy  is  then  distributed  to  the  other  energy  components  by  the 
pressure-velocity  interactions.  It  is  the  velocity-pressure  gradient  term  which  appears 
naturally  in  the  Reynolds-stress  transport  equations.  This  will  be  denoted  herein  as 

0,j  =  -  <UidpVdXj>  -  <u'j'dp'/3Xj>  (39) 

This  term  includes  both  the  pressure-strain  and  pressure-diffusion  terms.  Figure  12  gives 
the  velocity-pressure  gradient  term  for  all  six  stresses  in  the  M^.  =  0.2  and  M^.  =  2.5  cases. 
Because  of  the  minus  sign  in  Eq.  (36),  a  positive  value  in  Fig.  12  acts  to  reduce  positive  stress 
or  increase  a  negative  stress.  (For  consistency  of  notation,  let  (p-^2  =  (Puv 

The  terms  are  rather  similar  in  the  high  and  low  speed  results;  both  reduce 
<pu'’u’'>,  although  is  slightly  larger  for  the  low  speed  case.  Recall  from  Fig.  10  that 
<pu  "u  ">  is  also  similar  in  the  low  and  high  speed  flows.  The  terms  are  also  comparable, 
but  now  this  term  is  slightly  larger  in  the  high  speed  case.  However,  there  is  an  enormous 
’  hidden"  difference:  the  shear  stress  and  its  production  in  the  high  speed  case  is  only  a  small 
fraction  of  the  shear  stress  in  the  low  speed  case.  That  is,  the  ratios  -^y^,/<pu"v">  and 
increase  significantly  as  the  convective  Mach  number  rises.  The  "common 
wisdom  ”  in  turbulence  research  is  that  the  rate  at  which  the  shear  stress  is  reduced  by 
pressure-velocity  interactions  is  proportional  to  the  magnitude  of  the  shear  stress  and  its 
production  this  is  clearly  not  satisfied  across  the  range  of  convective  Mach  numbers 
considered  here. 

Another  large  difference  is  seen  in  0^.  This  term  is  responsible  in  part  for  energy 
transfer  from  <pu"u">  to  <pv"v”>.  The  rate  of  creation  of  <pv  "v">  by  0^^  in  the  high  speed 
layer  is  only  half  of  that  rate  in  the  low  speed  layer.  To  reiterate,  energy  is  extracted  from 
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the  mean  flow  as  <pu"u">  and  then  converted  to  <pv"v">  by  and  0^^;  <pv"v">  is 
required  to  generate  <pu"v">  via  and,  hence,  is  critical  to  the  turbulent  shear  stress.  The 
reduced  transfer  from  <pu"u">  to  <pv"v''>  disrupts  a  chain  of  events  which  normally 
produces  the  shear  stress. 

Are  these  statistical  results  consistent  with  the  swept  vortical  structures  of  turbulence? 
A  vortex  strongly  couples  the  fluctuating  velocity  components  normal  to  its  axis.  A 
spanwise  (z-oriented)  vortex  must  have  roughly  equal  magnitudes  of  u'  and  v'  irrespective 
of  w',  and  a  streamwise  (x-oriented)  vortex  has  roughly  equal  magnitudes  of  v'  and  w' 
regardless  of  u'.  This  is  consistent  with  the  turbulent  normal  stresses  <pv''v">  and  <p"w"w> 
in  Figure  10  which  are  roughly  equal.  Streamwise  fluctuations  u'  exist  in  this  flow  without 
being  strongly  coupled  to  v'  and  w'.  The  source  of  turbulence  energy  in  these  flows  is  the 
production  term  for  the  streamwise  fluctuations 

=  -  <pu"v">  dU/  dy  (40) 

The  streamwise  structures  are  unable  to  couple  this  source  of  turbulence  energy  to  the 
vertical  fluctuations  <pv"v">  which  are  required  to  generate  shear  stress.  Hence,  these 
statistical  details  of  turbulence  are  consistent  with  the  vortex  structure.  However,  the 
increased  suppression  of  <pu”v">  by  0^^,  cannot,  at  present,  be  linked  to  the  vortical 
structure. 

The  following  seem  to  be  significant  findings  from  the  analysis  of  simulation  data  at 
.M^  =  0.2  and  M,  =  2.5. 

1 .  In  the  Mj.  =  2.5  case,  the  velocity  profile  exhibits  a  sharp  transition  at  the  edge 
of  the  layer,  which  is  very  unlike  the  smooth  hyperbolic  tangent-like  profile 
in  the  low  speed  case. 

2.  The  turbulent  normal  stresses  are  very  anisotropic  at  M^.  =  2.5,  with  the 
streamwise  component  <pu"u">  much  larger  than  <pv”v  ">  and  <pw"w">. 

The  low  speed  turbulence  is  more  nearly  isotropic. 

3.  The  difference  between  the  full  stresses  <puju'j>  and  its  incompressible  part 
<p><u5uj>  is  small  when  there  is  no  density  stratification. 

4.  The  gross  compressibility  effect,  that  of  reducing  the  growth  shear  stress, 
may  be  explained  by  the  pressure-velocity  correlations  which  suppress  the 
shear  stress  and  which  do  not  rapidly  generate  <pv"v">  from  the  turbulence 
energy  in  <pu"u">.  Again,  <pv’’v”>  is  required  for  generation  of  the  shear 
stress  by  the  production  term  P^^. 

5.  The  swept  vortex  idealization  of  high  speed  turbulence  is  consistent  with 
these  statistical  results. 

The  statistical  analysis  is  incomplete.  The  velocity-pressure  gradient  terms  should  be 
decomposed  into  the  pressure-strain  and  pressure-diffusion  terms  to  gain  further  insight, 
and  then  the  pressure  strain  term  should  be  further  broken  into  the  so-called  return-to- 
isotropy  and  rapid  components,  if  possible.  It  may  also  be  important  to  determine  what 
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portion  of  the  reduction  in  shear  stress  at  high  Mach  number  can  be  attributed  to  specific 
statistical  mechanisms,  and  to  assess  the  possible  link  between  the  pressure-velocity 
statistics  and  the  swept  vortex  structure. 

The  Role  of  Dissipation  -  Proposals  have  been  put  forth  (e.g.,  Zeman^^)  concerning 
the  role  of  the  dissipation  rate  in  compressible  turbulence.  The  argument  is  that  dissipation 
rate  e  is  increased  due  *:o  shocklets  or  dilitation  which  reduces  the  turbulence  energy  k  and 
the  length  scale  L  -  k^'  e.  The  present  results  do  not  support  these  proposals. 

Firstly,  shocklets  are  virtually  non-existent  in  the  present  simulation.  Secondly,  the 
present  simulations  predict  the  reduced  spreading  rate  even  though  they  cannot  accurately 
predict  details  of  dissipation  phenomena.  The  sub-grid-scale  viscous  model  on  which 
dissipation  depends  is  an  eddy-viscosity  model  with  poor  accuracy,  and  the  accuracy  of  the 
differencing  scheme  becomes  increasingly  poor  at  the  high  wave  numbers  at  which 
dissipation  is  important.  Hence,  the  accuracy  for  dissipation  is  poor.  The  fact  that  these 
simulations  predict  the  decreased  mixing  at  high  despite  the  poor  accuracy  for 
dissipative  phenomena  strongly  suggests  that  dissipation  is  not  the  principal  cause  of  the 
reduced  shear  stress  in  high  Mach  number  turbulence. 

Limitations  -  Before  closing,  it  is  worth  reiterating  the  primary  assumptions  which 
may  compromise  the  accuracy  of  the  simulation  results.  In  most  cases  these  issues  are  not 
imposed  by  the  available  simulation  technology,  but  rather  by  the  limited  scope  of  a  Phase  I 
effort. 

1.  Most  simulations  are  of  a  temporally  evolving  shear  layer  in  which  the 
assumption  of  spatial  periodicity  is  made.  Streamwise  gradients  in  the  shear 
layer  thickness,  entropy,  and  the  evolution  of  the  dominant  structure,  which 
are  disregarded,  may  contribute  to  errors  in  the  results.  In  the  one 
comparison  between  similar  flows  that  were  computed  both  as  spatially  and 
temporally  evolving,  the  dominant  structures  were  very  similar.  Hence,  this 
is  not  felt  to  be  a  major  concern. 

2.  The  periodic  boundary  conditions  require  an  integer  number  of  the  dominant 
vortex  structures  to  fit  into  the  computational  domain.  This  means  that  there 
is  some  uncertainty  in  the  nature  of  these  structures,  which  depends  on  the 
size  of  the  computational  domain.  For  the  M^.  =  2.5  case  this  ambiguity  is 
small,  ±5®  in  sweep  angle,  and  does  not  materially  affect  the  conclusions  of 
this  study.  However,  there  may  have  been  a  significant  effect  in  the  case  of 
the  M^  =  4  simulation. 

3.  All  simulations  have  been  started  with  prescribed  initial  conditions,  and  the 
dependence  on  these  initial  conditions  is  large.  It  has  been  assumed  that  the 
random  initial  conditions  give  the  most  "natural "  well-developed  turbulence, 
but  this  is  not  proven. 

4.  These  are,  in  effect,  low  Revnolds-number  simulations  because  of  the  limited 
range  of  scales  which  are  resolved.  While  the  M^  =  0.2  solution  contained  a 
chaotic-looking  fine  scale  imposed  upon  the  dominant  vortical  structures,  the 
Mj-  =  2.5  case  was  more  nearly  dominated  by  the  large  scale  structure.  This 
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problem  may  have  biased  some  results,  especially  those  concerning  the  level 
of  anisotropy.  The  present  study  should  be  recognized  as  a  limited 
investigation  of  the  changes  in  turbulence  due  to  convective  Mach  number 
because  of  the  restricted  range  of  scales. 

5.  The  sub-grid-scale  stresses  are  treated  with  a  very  simple  model  which 
undoubtedly  has  poor  accuracy.  Hence,  the  small  scale  phenomena  related  to 
dissipation  will  be  represented  inaccurately. 

These  limitations  should  affect,  in  general,  details  in  the  simulation  results,  not  the  major 
conclusions.  Only  the  issue  regarding  the  form  of  the  initial  disturbance  can  have  a 
significant  effect  on  the  basic  findings  regarding  the  changes  in  turbulence  as  the  convective 
Mach  number  becomes  large. 

SUMMARY  AND  CONCLUSIONS 

The  objective  of  the  present  work  has  been  to  identify  major  changes  in  turbulence  as 
the  Mach  number  becomes  very  large.  The  present  effort  is  of  limited  scope  and  has  focused 
on  turbulence  in  an  ideal  gas,  although  it  is  recognized  that  many  factors  other  such  as  wall 
catalysis,  real  gas  and  rarefied  gas  effects  can  be  very  significant  in  hypersonic  flows.  The 
present  work  has  involved  theoretical  analyses  and  large  eddy  simulations  of  plane  free 
shear  layers 


The  present  results  provide  broad  support  for  the  idea  that  the  dominant  vortical 
structures  of  turbulence  become  swept  in  a  streamwise  direction  at  high  convective  Mach 
numbers.  One  analysis  is  widely  applicable  to  two-dimensional  shear  flows  in  the  limit  of 
the  Mach  number  tending  to  infinity.  It  indicates  that  entropy  constraints  restrict  vortices  to 
be  in  the  streamwise  direction,  which  minimizes  entropy  production.  Based  on  this  finding, 
a  simple  swept  vortex  model  for  planar  free  shear  layers  was  constructed.  This  model 
makes  predictions  about  the  sweep  angle  of  the  dominant  vortices  and  about  the  layer's 
growth  rate  over  the  full  Mach  number  range.  The  growth  rate  predictions  are  consistent 
with  experimental  results,  and  the  sweep  angle  predictions  are  generally  in  agreement  with 
simulation  results.  This  model  also  predicts  that  the  sweep  angle  is  large  enough  for  the 
flow  normal  to  the  vortices  to  remain  subsonic.  Hence,  shocks  will  generally  be  avoided  in 
plane  free  shear  flow. 

The  concept  behind  the  simple  vortex  model,  that  it  is  the  effective  velocity  normal  to 
the  vortical  structures  that  determines  the  shear  laveTs  growth,  is  not  unique  to  the  present 
effort.  However,  the  present  theorv  has  been  used  to  make  predictions  about  vortex  sweep 
angles  and  growth  rates,  which  are  confirmed  in  numerical  simulations  of  the  full  nonlinear 
vortex  development. 

The  underlving  mechanism  bv  which  high  speed  turbulence  avoids  shocks  is  not 
clearly  knowm.  The  present  result  suggests  that  it  may  be  a  straightforward  compressibility 
effect,  like  the  one  w’hich  effects  compressible  subsonic  spanwise  vortices.  Other  theories 
have  entertained  the  possibility  that  nascent  shocks  cause  spanwise  vortices  to  bifurcate  and 
sweep.  The  predictions  these  theories  make  about  the  shear  laver  growth  and  sweep  angle 
are  not  very  different.  However,  the  undedying  phvsics  are  different,  and  these  differences 
mav  become  critical  in  future  modeling  efforts  or  in  flows  other  than  the  plane  shear  layer. 
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Simulations  at  convective  Mach  numbers  of  M^.  =  0.2, 1 .48  and  2.5  were  evaluated.  The 
dependence  of  the  shear  layer's  growth  rate  upon  Mach  number  is  in  general  agreement 
with  experiment  and  the  vortex  model  theory.  The  structures  which  dominate  these  shear 
layers  are  spanwise  at  M,.  =  0.2  and  swept  about  75°  toward  the  streamwise  direction  for  the 
high  speed  flows.  The  effective  convective  Mach  number,  normal  to  the  structure,  is 
approximately  M^  -  0.6,  which  is  subcritical. 

Shock  waves  were  not  observed  in  these  randomly  initiated  simulations,  except  on  one 
very  brief  occasion.  The  absence  of  shocks  is  incompatible  with  theories  which  rely  on 
additional  shock-related  dissipation  to  account  for  the  reduced  mixing  in  high  speed 
turbulence. 

Significant  dependence  of  the  simulation  results  upon  the  perturbations  used  to  initiate 
the  calculations  was  observed.  The  flows  generated  bv  random  forcing  are  believed  to  be 
representative  of  "natural"  turbulence-initiating  disturbances.  Simulations  initiated  with 
specific  perturbation  modes  selected  from  stability  calculations  contained  structures  which 
were  spanwise  and  which  contained  strong  shocks.  These  results  are  probably  anomalous. 
However,  some  aircraft  may  preferentially  excite  specific  modes,  and  this  result  may  have 
important  implications.  This  result  clearly  indicates  that  the  starting  conditions  for 
simulations  must  be  treated  carefully. 

A  statistical  analysis  of  the  simulation  results  was  performed,  and  there  are  several 
significant  differences  between  the  M^.  =  0.2  and  2.5  cases.  It  can  be  argued  that  the  more 
important  among  these  are  related  to  the  change  of  sweep  angle  of  the  dominant  vortical 
structures.  For  example,  the  highly  swept  structures  cannot  rapidly  transfer  energy  of 
streamwise  fluctuations  into  normal  and  spanwise  fluctuations;  this  disrupts  a  critical  chain 
of  events  that  generates  shear  stress  (in  low  speed  flows). 

FUTURE  WORK 

There  are  at  least  three  areas  on  which  future  work  should  focus.  The  first  is  the 
resolution  of  certain  issues  raised  in  this  Phase  1  study.  Some  aspects  of  the  present  work 
were  based  on  unproven  assumptions  or  were  compromised  slightly  by  the  need  for 
expediencv.  Now  that  these  issues  are  reasonably  well  defined,  they  can  be  resolved.  One 
significant  issue  is  that  the  underlying  cause  of  vortex  sweep  is  not  known.  It  should  be 
possible  to  construct  careful  simulation  tests  to  identify  the  cause  of  this  phenomenon, 

A  second  aspect  of  any  future  work  is  the  consideration  of  more  complex  flows. 
Turbulent  flows  that  are  w’all-bounded,  affected  bv  mean  streamwise  vorticity,  three- 
dimensional  in  the  mean,  separating,  and  affected  by  shock  interactions,  for  example,  are  of 
great  importance  to  high  speed  aircraft.  Real  gas  and  rarefied  gas  effects  are  also  important. 
The  results  of  the  present  study  may  suggest  trends,  but  they  cannot  be  extrapolated  to 
these  other  flow  phenomena.  The  Mach  number  dependent  structural  changes  seen  in  the 
plane  shear  layer  may  or  mav  not  be  widely  significant  in  other  flows.  Even  simple 
phenomena  such  as  lateral  divergence  OW/dz  <  0,  in  the  present  notation)  ma\' 
significantlv  affect  the  characteristics  of  turbulence.  A  few  classes  of  flows,  which  embody 
critical  aspects  of  high  speed  turbulence  or  which  are  deemed  to  be  most  important 
technologicallv,  would  be  studied. 
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A  third  area  of  inportance  is  to  provide  for  improved  flow  prediction  methods.  Work 
in  the  first  two  areas  will  generate  understanding  of  the  underlying  flow  physics,  but  quite 
often  there  is  a  huge  gap  between  that  knowledge  and  the  statistical  models  required  for 
prediction  methods.  One  critical  aspect  of  this  work  would  be  to  identify  the  relationship 
between  turbulence  statistics  and  the  dominant  structure  of  turbulence,  assuming  such  a 
link  exists.  The  analyses  performed  in  the  present  study  should  provide  a  foundation  for 
developing  this  relationship.  The  present  work  suggests  that  the  vortical  structure  plays  a 
role  in  the  statistical  phenomena  which  must  be  modeled  in  existing  advanced  models.  If 
so,  then  some  means  of  solving  equations  for  the  structural  characteristics,  or  of  inferring 
structure  from  the  computed  turbulence  quantities,  must  be  developed.  The  logical 
conclusion  to  this  line  of  thinking  could  be  that  a  new  class  of  turbulence  models  is  required. 
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Figure  1.  Assumed  mixing  layer  growth  rate  versus  convtjtive  Mach  number  for  two- 
dimensional  vortex  structures,  ofaD^^c)- 


Figure  2.  Predicted  growth  rate  of  shear  layer  with  sweep  permitted,  t;3D(M.  ,0). 
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Figure  4. 
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Figure  3.  Sweep  angle  for  maximum  growth,  /3„  max^^c)- 


Temporal  growth  of  shear  layer  thickness  6/6q  as  function  of  non-dimensional 
time  tU/6o,  from  simulations. 
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Figure  5.  Equivalent  asymptotic  shear  layer  growth  rate  d6/dx,  from  simulations, 
normalized  by  d6/dx  at  M^,  =  0.2. 
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Figure  6(a).  Contours  of  side  view  of  instantaneous  spanwise  vorticity. 
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Figure  7(b). 


Figure  7(c). 


U  VELOCITY 


Contour  plots  for  =  2.5  case;  side  view  of  stream  wise  velocity,  showing 
streamwise  inclination. 


Contour  plots  for  M^.  =  2.5  case;  oblique  view  of  streamwise  vorticity  from 
above  shear  layer. 
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NORMALIZED  PRESSURE 
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Figure  10(b). 


Figure  11(a). 


Turbulent  stresses  <fiu'-u'l>  for  2.5  case.  (Notation  is  a=  <pu”u">, 
+  =  <pv"v”>;  0  =  <pw"w">;  o  =  <pu''v">;  A  =  <pu’'w">;  x  =  <pv"w">.) 
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Decomposition  of  the  total  turbulent  stresses  (T)  into  incompressible  (1)  and 
purely  compressible  (C)  parts,  defined  by  Eq.  36  and  37,  for  <pu"u''>  at 
M  =0.2. 


STRESS  COMPONENTS 


Figure  11  (b). 


Figure  11(c). 


Decomposition  of  the  total  turbulent  stresses  (T)  into  incompressible  (I)  and 
purely  compressible  (C)  parts,  defined  by  Eq.  36  and  37,  for  <pu”u">  at 
M,  =  2.5. 
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Decomposition  of  the  total  turbulent  stresses  (T)  into  incompressible  (I)  and 
purely  compressible  (C)  parts,  defined  by  Eq.  36  and  37,  for  <pu”v”>  at 


